DO-TH-02/22 



Nonequilibrium evolution of $ 4 theory in 1 + 1 dimensions in the 2PPI formalism 
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We consider the out-of-equilibrium evolution of a classical condensate field and its quantum fluc- 
tuations for a $ 4 model in f + 1 dimensions with a symmetric and a double well potential. We 
use the 2PPf formalism and go beyond the Hartree approximation by including the sunset term, 
fn addition to the mean field cj>(t) = (<3>) the 2PPf formalism uses as variational parameter a time 
dependent mass M 2 (t) which contains all local insertions into the Green function. We compare our 
results to those obtained in the Hartree approximation, fn the symmetric $ 4 theory we observe that 
the mean field shows a stronger dissipation than the one found in the Hartree approximation. The 
dissipation is roughly exponential in an intermediate time region. In the theory with spontaneous 
symmetry breaking, i.e., with a double well potential, the field amplitude tends to zero, i.e., to the 
symmetric configuration. This is expected on general grounds: in I + I dimensional quantum field 
theory there is no spontaneous symmetry breaking for T > 0, and so there should be none at finite 
energy density (microcanonical ensemble), either. Within the time range of our simulations the 
momentum spectra do not thermalize and display parametric resonance bands. 

PACS numbers: 03.65.Sq, 05.70.Fh, 11.30.Qc 
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I. INTRODUCTION 



There has recently been considerable activity in inves- 
tigating the nonequilibrium evolution of quantum field 
theory beyond the large- N approximation. In particular 
there has been formulated a systematic approach 

(2PI-NLO) in which all next-to leading order contribu- 
tions in the 1/N expansion are included, using the CJT 
or 2PI formalism |3( generalized to nonequilibrium evo- 
lution in Ref. £|, using the closed time path (CTP) or 
Schwinger-Keldysh 5] formalism. A similar approxima- 
tion including some but not all contributions of next-to- 
next-to leading order is the bare vertex approximation 
(BVA) 6]. Numerical simulations have been performed 
mostly in 1 + 1 dimensions, for the classical or quantum 
<I> 4 theory with a symmetric Q, 0, H, 0, and with a 
double well ^lj potential. A first simulation in an O(N) 
model for N ^ 1 in 3 + 1 dimensions has just appeared 
[l2|. for a symmetric potential. The analysis of models 
with N 7^ 1 and with spontaneous symmetry breaking, 
i.e., a Mexican hat potential, should be very important 
for understanding the role of the Goldstone modes and 
their influence on the phase structure of the theory in a 
certain approximation. 

A more modest step beyond leading order large TV has 
been taken in Ref. [13j , where the Hartree approximation 
was used in an O(N) model in 3 + 1 dimensions. This 
approach includes some terms of the nonleading orders, 
but not all of them. Our investigation made evident the 
role of parametric resonance in the system of "sigma" 
and "pion" modes, and the role of the Goldstone modes 
in stabilizing the evolution in the regions where the equi- 



* Electronic address: 
t Electronic address: 



baacke@physik.uni-dortmund.de 



andreas.heinen@uni-dortmund.de 



librium effective potential is complex, features still to be 
investigated in the systematic 1/N motivated approxi- 
mations. 

An approach for including nonleading orders in a self- 
consistent resummation scheme has been proposed some 
time ago by Verschelde and collaborators |14lllfi| . the so- 
called 2PPI expansion. Like in the 2PI or CJT formalism 
the mean field and the internal Green functions are de- 
termined self-consistently. As in the 2PI formalism the 
equations of motion follow from an effective action, which 
here is a functional of a mean field (f> = (<&) and an effec- 
tive mass M. . In contrast to the 2PI approach only local 
insertions into the Green function are resummed, so the 
Green function is, in all orders, a functional of the local 
mass term M. that in general will depend oni= (i,x). 
In particular this Green function is different from the 
physical Green function. As far as the resummation is 
concerned the approach is less powerful: one has to in- 
clude more diagrams if one wants to reach the same order 
in a loop or 1/N expansion as in the 2PI expansion. The 
fact that the propagators have a simpler structure may be 
a disadvantage as by this structure the approximation is 
less flexible than 2PI. On the other hand the calculations 
are technically less involved, in particular the formalism 
does not require ladder resummations which complicate 
the renormalization of the 2PI approach 0] . In one-loop 
order the approach is equivalent to the Hartree approxi- 
mation. Recent progress in the 2PPI formalism includes 
the demonstration or renormalizability [171 [18J and some 
finite-temperature two-loop calculations in 3 + 1 dimen- 
sions. An interesting result was that for N = 1 and 
for N ^ 1 |2(J the order of the phase transition between 
the spontaneously broken and symmetric phases becomes 
second order in the 2-loop approximation, while it is first 
order in the Hartree approximation. The results for the 
2PPI expansion have been compared to exact results in 
[2l| for the anharmonic oszillator; even more recently j2^ 
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the two-loop approximation has been compared, in 1 + 1 
dimensions, with exact results of the Gross-Neveu model. 

In this paper we will present the formulation and 
some numerical results for the noncquilibrium evolution 
of the mean field and the self-consistent mass in the 2PPI 
scheme, applied to $ 4 theory in 1 + 1 dimensions. We go 
beyond the one-loop (Hartree) approximation by includ- 
ing the sunset graph, which represents the full two-loop ( a ) 
contribution in this formalism. We explicitly formulate 
a conserved energy functional which is used to monitor 
the numerical accuracy. As this is the first investigation 
of the 2PPI formalism at two loops out of equilibrium 
we do not attempt a detailed study; rather we aim at 
presenting the main new features of this approach, as 
well as compared to the one-loop Hartree approximation 
and as compared to the other approaches (2PI-NLO and where 
BVA) mentioned above. We consider both the case of the 
symmetric $ 4 potential and the case of the double well 
potential which displays spontaneous symmetry breaking 
on the classical as well as on the one-loop level. 

The plan of the paper is as follows: In section 2 we 
formulate the model and present the 2PPI formalism as 
applied to the system out of equilibrium. In section 3 
we specify the two-loop approximation by giving the ex- 
plicit expressions for the basic graphs, the equations of 
motion, the conserved energy and by discussing the ini- 
tial conditions and renormalization. In section 4 we give 
details of the numerical implementation. The numerical 
results are presented and discussed in section 5. We end 
with conclusions and an outlook in section 6. The paper 
is completed by two appendices. 
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FIG. 1: Examples for two particle point reducible (2PPR) 
and irreducible (2PPI) diagrams: (a) a diagram which is 2PR 
and 2PPR (b) a diagram which is 2PR but 2PPI; solid lines: 
internal propagators; dashed lines: external fields 4>. 



M 2 {x) 



V(x) + ^A(x). 



(2.3) 
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So in contrast to the 2PI formalism this Green function 
is not a variational object, it is a functional of A4 2 (x). It 
is not the physical Green function. 

The 2PPI formalism in its original form is based on 
the action 



r[0,A] = 5 class [0] + r 2PP1 [</>,M 



II. FORMULATION OF THE MODEL 

We consider the $ quantum field theory defined by 
the Lagrange density 

C = -d B $a"$--m 2 $ 2 -- $ 4 . (2.1) 
2 M 2 24 y J 

If m 2 > we refer to as the symmetric theory, and with 
m 2 < which we refer to as theory with spontaneous 
symmetry breaking. These terms relate to the classical 
theory and do not imply the occurrence of spontaneous 
symmetry breaking in the quantum field theory. 

The 2PPI formalism proposed by Verschelde and Cop- 
pens [l4l fill is based on an effective action which is for- 
mulated in terms of a mean field <f> and a local insertion 
A. It is the Legendre transformation of a generating 
functional with a source J(x) for the field 4>(x) and an- 
other local source K(x) for 4> 2 {x). Here lies the difference 
to the well-known CJT formalism, where one introduces T[4>,M.' i 
a bilocal source K(x,x') for a Green function G{x,x'). 
Graphically the 2PPI scheme resums all local insertion 
into a Green function which in all orders remains a gen- 
eralized free particle Green function, or Green function 
in an external field, i.e. 



(2.4) 

Here P 2PPI [</>, M 2 } is the sum of all 2PPI graphs; these 
are defined as graphs which do not decay into two parts of 
two lines joining at a point are cut. In the 2PI formalism 
one includes into the analogous T 2PI all graphs which do 
not decay into two parts if any two lines are cut. In or- 
der to visualize the difference we show examples of 2PR 
and 2PPR graphs in Fig. ^ It should be emphasized 
that in the 2PPI formalism the lines in the graphs are 
Green functions with the variational mass term A4 2 (x) 
as defined via Eq. I|2.2|) . while in the 2PI formalism the 
internal lines refer to the variational Green functions in- 
duced by the bilocal sources. When comparing the sets 
of irreducible graphs in both formalisms one has to take 
into account this difference in the meaning of internal 
lines. 

The insertion A(x) is given by 



A(x) 



ST 



2PPIr 



,M 2 



8M 2 {x) 



(2.5) 



As we have stated before it is simpler to formulate the 
action in terms of 4> and M 2 . We solve Eq. (|2.3(l with 
respect to A and insert this into Eq. 1)2. 4f> . Using the 
explicit form of S c ia SS we obtain 



d 2 x 



^d^(x)d^(x)~^M 2 (x)^ 2 (x) 



G~\x, x') = i [□ + M 2 {x)} 8(x - a;') , 



(2.2) 



+ ^Jd 2 x[M 2 (x)-m 2 ] 2 
+T 2PPl \^M 2 ] . 



(2.6) 
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FIG. 3: Bubble and sunset diagrams. We display the lead- 
ing diagrams in the 2PPI action T 2PPI : (a) the bubble dia- 
gram; (b) the sunset diagram; solid lines: internal propaga- 
tors; dashed lines: external fields (f>. 



FIG. 2: The closed time path in the complex t plane. 



One easily checks that the equations of motion obtained 
by varying this action with respect to 4>(x) and A4 2 (x) 
take the form 







□0 + M 2 (x)<j)(x) - ^(t) 3 (x) 



5T 



2PPIr 



A 



5<t>{x) 
5T 2PPl \6,M 2 ] 



M 2 (x)= m 2 + -<p 2 (x)-X 5M2{x) 



(2.7) 
(2.8) 



The latter equation being identical to Eq. 1)2. 3[) . we will 
refer to it as gap equation. 

Here we will consider states of the system that are 
spatially homogeneous; so in Eqs. Ij2.7|l and (|2.8(l the 
arguments x = (t, x) should be replaced simply by t. 
Furthermore, in Eq. (|2.6|l the space integration simply 
gives a volume (length) factor, and in the nonequilibrium 
formalism the time integration should be replaced by the 
closed time path (CTP) displayed in Fig. 

The equation for the Green function separates into 
space and time dependence. Using the homogeneity of 
the state in space the Green function can be written as 



G(t, t'; x, x) 



^e zp{x ~ x,) G{t,t'; P ) . (2.9) 
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As the equation is local in time the Green function 
G(t,t';p) can be expressed in terms of mode functions 



A 



A 



m 2 = M 2 (0) =m 2 + -«^(0) + -A(0) 



This means that /(0,p) = 1 and /(0,p) 



(2.12) 



-iujp with 

Lo p = ^Jp 2 + tjiq. This defines an initial Fock space. We 
continue the discussion on initial conditions in section 

In the CTP formalism § one uses Green functions 
with different time orderings. The Green function 
G(t 7 t';p) as defined above is identical to the Green 
function G>(t, t';p) with normal time ordering; the 
anti-time-ordered Green function G<(t, t';p) is given by 
G < (t,t';p) — G>(t',i;p). In the explicit formulae one 
can use this identity in order to express all relevant Green 
functions by G(t,t',p). 



III. ONE-LOOP AND TWO-LOOP 
CONTRIBUTIONS 



Having established the general formalism we can now 
discuss the leading terms in a loop expansion. The two 
relevant graphs are depicted in Fig. |3 graphs a and b. 
The leading bubble diagram is the "log det" contribution. 
It leads to the tadpole insertion into the Green function. 
The next diagram represents the only contribution on 
the two-loop level. We will separately discuss these two 
contributions in the following. 



G(t,t';p) 



2cj, 



[f(t,p)r(t', P )Q(t-t') 

+f(t',p)r(t,p)Q(t'-t)} (2.10) 



where lo p is defined below, and where the mode functions 
f(t,p) satisfy 



f(t,p)+[p 2 +M 2 (t)] f(t, P ) = 0. 



(2.11) 



Here we choose the initial conditions for the mode func- 
tions at t = as for wave functions of free particles with 



A. The bubble diagram 

The bubble diagram defines the leading Hartree con- 
tribution. It is independent of (f> and so does not yield an 
explicit contribution to the equation of motion for <fi. Of 
course it enters, indirectly, via its contribution to M. 2 . 
Its contribution to T 2PFl ((j), M 2 ) is given by 



T^[^M 2 } - ^Trln [G^iM 2 ]] 



(3.1) 
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Its functional derivative with respect to M 2 is given by naive quantum energy defined by 



6M 2 {t) 



dp 



I'kIlo, 



\f(t, P W 



(3.2) 



The energy density can be derived [23j by considering 
a variation of the action under t — > t + r(t), which in- 
duces S(f)(t) = <i>(t)T(t), 5<j>{t) = 4>(t)r(t) +<p(t)i(t) and 
SM 2 (t) = T(t)dM 2 {t)/dt. The one-loop action only de- 
pends on A4 2 (t). One then finds that the contribution of 
the bubble graph to the energy is defined by the relation 



d£«(t) 



ST^[(I),M 2 } dM 7 



dt 



5M 2 {t) 



dt 



-G(t,t;p) — 



(3.3) 



This equation can be integrated explicitly; indeed one 
checks easily, using the mode equation (|2.11(l . that the 



dp 



2tt2w 



+ [ P 2 +M 2 (t)]\f(t 7 p)\ 2 . (3.4) 



is consistent with the defining equation 1)3.3(1 . This is 
of course well-known. If only this one-loop contribution 
is included, the approximation is referred to as Hartree 
approximation. 



B. The sunset diagram 

The unique two-loop contribution to r 2PPI is the sun- 
set diagram which in the CTP formalism is explicitly 
given by 



T {2) [(j),M 2 }=i— I dxdx' I dt<p{t) \ dt' G%{t,t';x,x')<j){t') 



(3.5) 



where the t and if integrations are over a CTP contour and where Gp is the path-ordered Green function. 
The functional derivative with respect to <j>(t) is given by 



-S(t) 



with 



5 ( *) = -4/ *Wn(frW£i 

JO. J o — i \ / \ t>—-\ 



l[G(t,t';pc)-l[G(t',t,p e ) 



(3.6) 



(3.7) 



This contribution to the equation of motion for <f>(t), an amputated sunset diagram is represented graphically in 
Fig. EH. 

The functional derivative with respect to A4 2 (t) is given by 



6M 2 (t) 



1 



2 

2 



x [G(t,t';p 3 )-G(t',t;p 3 )} 

x [G(t',t , '; Pl )G(t'X-, P2 )G(t,t' , ;P3)-G(t , ',t , ; Pl )G(t , ',t'; P2 )G(t' , ,t-, P3 )} 



(3.8) 



r 



This graph which contributes to the gap equation is de- Considering again a variation t — ► t + r(t) one finds 

picted in Fig. EJs. It is a tadpole diagram with fish the contribution of the sunset term to the energy to be 
insertion. 
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FIG. 4: Two-loop graphs in the equations of motion, (a) 
the amputated sunset diagram which appears in the equation 
of motion l|2.7|l for the mean field <f>; (b) the tadpole diagram 
with fish insertion which contributes to the gap equation 12.81 1 ; 
in both diagrams the solid dots indicate the external time, 
the time variables of the other vertices appear in the internal 
integrations; solid lines: internal propagators; dashed lines: 
external fields d>. 



defined by 

dE^ 2 \t) 
dt 



= <j>(t)S(t) + l d -^A^(t). (3.9) 



As far as we see this expression cannot be integrated 
explicitly; but the relation can be integrated numerically 
to obtain E {2 \ 



E^ 2 \t) = / dt 1 
Jo 



(3.10) 



D. Renormalization and initial conditions 

There is a wide choice of initial conditions for the sys- 
tem. So one may choose an initial mean field 0(0) and 
one may modify the Green function by including con- 
tributions from the kernel of G~ 1 {x,x l ). In the one- 
loop approximation the latter possibility is equivalent to 
choosing initial ensembles for which the modes f{t,p) are 
populated, or to Bogoliubov rotations of the initial Fock 
space. In the two-loop approximation this simple particle 
picture is no longer appropriate. 

However the choice of initial conditions is not en- 
tirely arbitrary because of initial singularities EH Up- 
starting with some nonzero value of 0(0) and with tuq a 
solution of the gap equation no initial time singularities 
are encountered. So such a choice is a physically accept- 
able one. In order to solve the gap equation at t = 
we have to know the contributions AW(0) and A( 2 )(0). 
A^ 2 ) has already be defined such that it vanishes at t = 0; 
by this choice we erase the memory of the past. AW is 
given by an integral over the fluctuations, so it does not 
vanish. Furthermore it is divergent and so we have to 
discuss renormalization. 

In <I> 4 theory in 1 + 1 dimensions there is only one prim- 
itive divergence, the one of the tadpole graph. Renormal- 
ization reduces, therefore, to making a shift in the tad- 
pole term which can be absorbed by a shift in the mass. 
Using dimensional regularization we rewrite the tadpole 
contribution A^, see Eq. (|3.2|l . as 



C. Equations of motion and energy 

Having determined the functional derivatives of the ac- 
tion with respect to cj>(t) and Ai 2 (t) we can explicitly 
write down the equations of motion in the two-loop ap- 
proximation: 







$(t) + M 2 (tMt)-±<p(t) + S(t) (3.11) 



M 2 (t) = m 2 + ~ \(f) 2 (t)+A^(t)+A^{t)\ (3.12) 

We define the "classical" energy as the zero-loop expres- 
sion 

2*<°>(t) = i^(t) + i^W(*)-^(*) 



-±-[M 2 (t)-m 2 ] 2 



We have defined the one-loop and two-loop contributions 
in Eqs. (|3.4|l and (|3.10|) . One can check, using the equa- 
tions of motion, that the total energy 



E (0) , E (l) , E {2) 



(3.14) 



is conserved. 



A«(*)=A«(*) + -M^- 7 + ln 



47T/i 2 



Air [e 

with the finite part of A^ defined as 
f d P 



2tt2w 



v 



1/(^)1-1 



(3.15) 



(3.16) 



This is the expression used in the numerical computation. 

Including a mass counter term the gap equation now 
takes the form 



M 2 (t) = m 2 + Sm 
A 
+ 2 



2 (t) + A«(i) 



(3.13) Choosing 



5m 



1 f 2 , Ann 2 

-— --7 + m— V 
47T e mi 



A f 2 , iTTfl 2 

7 + m 



(3.17) 
A( 2 >(i) . 

(3.18) 



8ir |_ e ;// 
the finite gap equation takes the form 

M\t) =m 2 + Sm 2 n + \ \tf(t) + A«(i) + A( 2 )(t) 



(3.19) 



6 



with 



DISCUSSION OF THE RESULTS 



A TO 2 

87r mi 



(3.20) 



Initial conditions and renormalization are equivalent to 
those in the one-loop approximation, which facilitates a 
comparison between the two-loop 2PPI and the Hartree 
approximation. 



IV. NUMERICAL IMPLEMENTATION 

In the 2PPI approximation the Green function factor- 
izes and so one can work with mode functions. This 
considerably facilitates the numerical computation of the 
"memory" integrals introduced by the sunset graph. In 
particular one has to store only functions of one time 
argument, of course still for all times and all momenta. 
The storage requirements grow only linearly with time, 
so the evolution can be followed for relatively long times. 
Furthermore the differential equations are ordinary dif- 
ferential equations that can be solved precisely using a 
Runge-Kutta algorithm. This can be important if one 
has to trace parametric resonance phenomena. Of course 
if the approximation itself is poor these numerical advan- 
tages are useless. Still, the possibility of doing the calcu- 
lations with good precision allows to study the quality of 
the approximation reliably, including its possible short- 
comings. 

The time integration was done in steps of At = 0.001 
to 0.005. The Wronskians of the mode functions were 
constant with a relative precision of 10 -8 . For the mo- 
mentum cutoff, which is a cutoff of a convergent integral, 
we have chosen p max = 20. As one can see from the 
momentum spectra, this is a rather generous choice. It 
should be mentioned that momentum conservation leads 
to momenta that can be beyond the cutoff. This is a 
problem that can hardly be avoided, and a relatively 
large momentum range should make such "losses" tol- 
erable. A more serious problem is the momentum grid. 
We observe parametric resonance [35| , and this leads to 
amplitudes that vary strongly in time and momentum. 
In the typical large- A studies this fact has lead the vari- 
ous groups to choose much finer grids with several thou- 
sand momenta. This is not possible here, we think that 
the essential features of the low momentum region with 
parametric resonance and/or exponential growth subsist 
with a less refined grids. This concerns in particular the 
self-stabilization of the system in the classically unstable 
regions. We have chosen Ap = 0.05, i.e., a grid of 400 
equidistant momenta. In principle such a grid can lead 
|27| to "lattice artefacts" , corresponding here to a lattice 
size L = 2-n j Ap = 407r in inverse mass units. Indeed we 
do not observe any phenomena that suggest such arte- 
facts. The choice of Ap is also discussed in Appendix 
B. 



We have performed several simulations for the case of 
a symmetric $ 4 potential and for a double well potential 
which classically leads to spontaneous symmetry break- 
ing. The initial configuration has been, in all cases, a 
mean field 4> different from its classical expectation value, 
and a quantum ensemble corresponding to the ground 
state of a Fock space characterized by an initial mass 
mo = A^(0). We have obtained results for the time evo- 
lution of the mean field 4>(t), for the self-consistent mass 
M. 2 (t) and for the energy. The relative importance of the 
two-loop contributions can be seen in their contribution 
to M 2 . In all cases we have compared the evolution with 
the one obtained in the Hartree approximation. 



A. Results for the symmetric <E> 4 potential 

In Figs.[S]and|S]we display our numerical results for the 
time evolution of the mean field (Figs. and[|Ji), of the 
dynamical mass M 2 (t) (Figs.[5|D andEb), of the sunset 
contribution S(t) in the classical equation of motion l|2.7|) 
(Figs.EJ; and|SJ:) and of the classical and quantum parts 
of the energy (|3.14|) (Figs. [SJi andHJl). In the latter 
diagrams we define the classical energy as the standard 
expression 

Ecl = 1 -4> 2 + 1 -(m 2 + Sm 2 a )cb 2 + ^. (5.1) 

Indeed the repartition between classical and quantum en- 
ergy is to some extent arbitrary in a self-consistent frame- 
work where, e.g., M 2 contains classical as well as quan- 
tum parts. In Figs.|SJi,b and[^L,b we also display the time 
evolution in the one-loop or Hartree approximation. 

We observe the following characteristic features: after 
an initial period of time in which the field amplitude stays 
roughly constant and close to the Hartree time evolution 
a period of effective dissipation sets in. For small initial 
amplitudes evidently the dissipative phase ends and the 
mean field reaches a roughly constant amplitude of oscil- 
lation, again. For large initial amplitudes such a "shut 
off" is less evident. A closer investigation shows that ini- 
tially the quantum modes build up until the sunset dia- 
gram becomes important. From then on the Hartree and 
two-loop evolutions differ substantially. The increase of 
the sunset diagram triggers dissipation, until the sunset 
diagram again becomes small due to the decrease of the 
external fields. Once the sunset diagram has lost its im- 
portance the amplitude of oscillation of the classical field 
becomes roughly constant again. This is seen in partic- 
ular in Fig. |Sp, where, due to a relatively small initial 
amplitude, the quantum modes and therefore the sunset 
diagram are less important than for large initial ampli- 
tudes (or energy densities), as, e.g., in Fig.[5Ji. We have 
not followed the evolution at really large times. So we 
cannot decide between a constant and a slowly decreasing 
amplitude as found, e.g., in the large- A case |3l| . 
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FIG. 5: Time evolution for the symmetric <E> 4 potential; Parameters: m 2 = 1 (symmetric potential), A = 6, 0(0) = 1.2; (a) 
evolution of the mean field; (b) evolution of the effective mass M 2 ; (c) evolution of the sunset contribution to the equation 
of motion of 0, 12.71 : (d) evolution of the energy; in (a-c) the solid lines relate to the the two- loop 2PPI approximation, the 
dashed lines to the one-loop or Hartree approximation; in (d) the dashed line is the classical energy 15. 1| . the dotted line is the 
quantum energy. 



The total energy, displayed in Figs.EJi andHJl, is con- 
stant as it should. Numerically this is the case within 
five significant digits or better; here was obtained 
by Runge-Kutta integration of Eq. t\'d. 9fl . 

We also present, in Fig. [7| typical momentum spec- 
tra. We have chosen the simulation with 0(0) = 1.2 and 
show the spectra for an early time t = 10 and at the 
end of the simulation. Along with the results for the 
two-loop approximation we display those for the Hartree 
approximation. Obviously the spectrum evolves more 
strongly for the two-loop approximation. At late times 
it shows the typical features of a parametric resonance 
band [H |H H<| Hj . Above this band the spectra drop 
to small values and decrease to zero. It should be ev- 
ident that our momentum cutoff of p max = 20 will be 
sufficiently high even for the multiple integrals. On the 
other hand the numerical integration cannot take into 
account the finer details of the spectra, in particular at 
later times. Of course to some extent these details are 
washed out if averaged over time. Still, to some extent 
the finer details of the time evolution of <p(i) show some 
dependence on the choice of Ap. However, neither here 
nor in the case of the double well potential the qualitative 
features are affected by these details. 

We have restricted our presentation to one single cou- 



pling parameter A = 6. We have performed simulations 
for smaller values of A, as well; for such values of A the 
time evolution is stretched; the dissipation sets in later 
and extends over a larger span of time. For A = 1 the 
dispersive phase extends to typically t = 300. The gen- 
eral, qualitative, characteristics of the time evolution are 
similar. 



B. Results for the double well potential 

The numerical simulations for the double well potential 
are presented in Figs. 151 Itll and ITU1 We take the coupling 
A = 1 and m 2 = —1/6, so that the classical minimum of 
the double well potential is at — v — 1- We consider 
initial values 0(0) equal to 1.5,1.4 and 1.2. Classically 
the system can cross the barrier between the two minima 
for 0(0) > y/2. The first of our initial values is above this 
critical value, the second one is slightly below it. In the 
Hartree approximation the system evolves as expected 
from this classical consideration. In the two-loop 2PPI 
approximation the system evolves towards the symmet- 
ric phase where the system oscillates around = in 
the later stages of evolution. The transition between a 
motion in the region of the classical minimum and = 
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FIG. 6: Same as Fig.EJfor 0(0) = 0.6. 
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FIG. 7: Momentum spectrum at time t = 10 and t = 50 for 
the parameter set of Fig. 



is accompanied by an increase of the sunset contribu- 
tion and of the mass A4 2 . The transition happens early 
for (f)(0) — 1.4. In this case we are just below the critical 
value, one sees that the transition towards <f> ~ happens 
at a time where the sunset contribution is still small and 
where M. 2 only slightly deviates from its Hartree value. 
For <j}(0) = 1.2 we are deeply in the well. Here it takes 
a long time before the evolution towards (f> ~ sets in. 
If we start with values (f)(0) even nearer to the classical 
minimum = 1 the transitions happens at even later 
times and we expect the discretization of the momen- 
tum spectrum to affect our results so as to make them 
unreliable. 

In Fig. ^2 we display momentum spectra for \f(t,p)\ 2 — 
1 for the simulation with <f>(0) = 1.4 at t = 20 and at 
t = 50, along with the spectra obtained in the Hartree 
approximation. For the two-loop simulation the spec- 
trum at t = 20 is characterized by a strong peak at low 
momentum, which apparently is due to a passing of M 2 
to slightly negative values. At t = 50 the effective mass 
of the modes is positive, the spectrum shows a charac- 
teristic band as typical for parametric resonance. 

In all simulations the Hartree approximation displays a 
rather clean periodicity which signals a strong coherence 
between the evolutions of the classical field and of the 
quantum modes. This effect is much stronger than in 
3+1 dimensions, in the Hartree ^| or large- TV [30l l32]| 
approximations. 
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FIG. 9: Same as Fig. for 0(0) = 1.4. 
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FIG. 10: Same as Fig.|Ufor 0(0) = 1.2. 
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FIG. 11: Momentum spectrum at time t = 20 and t = 50 for 
the parameter set of Fig. 



The nonequilibrium evolution of <I> 4 quantum held the- 
ory with a double well potential has been studied recently 
by Cooper et al. [Ill], with somewhat different initial 
conditions. These authors hnd a transition towards a 
symmetric phase in the 2PI formalism extended to next- 
to leading order in 1/7V (2PI-NLO), while in the bare 
vertex approximation (BVA) the system remains in the 
broken phase for initial energy densities below some crit- 
ical value. The exact theory has no phase transition at 
finite temperature and, therefore, is not expected to have 
one at finite energy density |33|. 



VI. SUMMARY, CONCLUSIONS AND 
OUTLOOK 

In this paper we have considered the out-of-equilibrium 
evolution of a classical condensate field 4> — ($) and its 
quantum fluctuations for a <I> 4 model in 1 + 1 dimensions, 
with a symmetric and a double well potential. Our 
investigation was based on the 2PPI formalism in the 
two-loop approximation. We have generalized the 2PPI 
formalism to nonequilibrium quantum field theory. In 
order to find the main features of this approximation we 
have performed a first set of numerical simulations and 
compared the results to the ones obtained in the Hartree 
approximation. 



We summarize our results as follows: 



11 



In the symmetric $ 4 theory we observe that the mean 
field shows a stronger dissipation than the one found in 
the Hartree approximation. The dissipation is roughly 
exponential in an intermediate time region. This dissi- 
pation is obviously related to the sunset contributions. 
As these involve the mean field amplitude they become 
unimportant when the amplitude goes to zero. There- 
fore, for later times the system seems to develop a stage 
of weak dissipation. However, we have not extended our 
study to "late" times in the sense of an asymptotic anal- 
ysis of the evolution. 

In the theory with spontaneous symmetry breaking, 
i.e., with a double well potential, the field amplitude 
tends to zero, i.e., to the symmetric configuration. This 
is expected on general grounds: in 1 + 1 dimensional 
quantum field theory there is no spontaneous symmetry 
breaking for T > 0, and so there should be none at finite 
energy density (microcanonical ensemble), either. 

We observe in both cases that parametric resonance 
phenomena are important, and that the momentum spec- 
tra show no sign of thermalization. In contrast to the 2PI 
approximation the interaction between the modes is via 
the spatially homogenous (zero momentum) mass term; 
so there is no direct momentum exchange between the 
modes via a Schwinger-Dyson equation and the analysis 
of Ref. [34| concerning thermalization in the 2PI for- 
malism therefore does not apply. Our numerical analysis 
does not allow definite conclusions about thermalization 
at later times. 

In conclusion we have shown that the 2PPI formal- 
ism can be generalized to nonequilibrium quantum field 
theory and that the simulations in the two-loop approx- 
imation in 1 + 1 dimensions show sizeable differences 
when compared to the Hartree approximation. Both the 
stronger dissipation and the correct symmetry structure 
overcome obvious deficits of the Hartree approximation. 
We therefore think that it is worthwhile to further inves- 
tigate the properties of this approximation, in and out of 
equilibrium. 

Obvious generalizations of this investigation include 
the analysis of an O(N) model with N > 1 in 1 + 1 di- 
mensions and analogous studies in 3 + 1 dimensions. The 
technical requirements for such simulations are consider- 
ably reduced when compared to the 2PI formalism in the 
analogous approximation, due to the factorization of the 
Green functions; moreover the problem of renormaliza- 
tion in 3 + 1 dimensions has been solved in equilibrium 
[l7tll8j . In the 2PI approach the three- loop renormaliza- 
tion has been considered in jlfj for the mean field = 
case; an analysis of renormalization beyond the Hartree 
approximation is still lacking for <f> ^ 0. 

We feel that it is very important to accompany the nu- 
merical simulations of nonequilibrium systems in various 
formalisms and approximations by equivalent analyses 
for systems in thermal equilibrium. Such analyses are 
still lacking entirely. 
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APPENDIX A: COMPARISON BETWEEN 2PI 
AND 2PPI 

In this section we are giving some comments on the 
differences between the 2PI (two particle irreducible) and 
the 2PPI (two particle point irreducible) formalism at the 
two-loop level. 

The 2PI effective action reads 

T[(f>,G\ = S[<p] + iiTrlnG -1 + i*T5r 

+r 2 [4>,G\, (A.l) 

where iD~ x (x, x') = -(D+m 2 )6(x-x')-^(f) 2 (x)S(x-x') 
is the classical propagator. If we truncate Y2 including 
two-loop order terms we have 



-.(2) r 



i,G] = 




(A.2) 



(2) 

A variation of Y 2 ' with respect to G gives the self energy 
E^(x,x') = 2/- 



: 5T™[<j>,G] 



SG{x,x') 




(A.3) 



The two point function G fulfills the Schwinger-Dyson 
equation Q 

iG-\x,x') =iD- 1 (x,x')-ii: (2) {x,x') (A.4) 

and is a variational parameter of the formalism. 

The relevant formula for the 2PPI formalism in the 
two- loop approximation are given in section Till Bl Note 
that although the sunset contribution to the effective ac- 
tion in both the 2PI and 2PPI approach is depicted by 
the same diagram, they have different implications. We 
think that it is instructive to compare some implicit 1PI 
graphs of both approximations to emphasize the differ- 
ences. These 1PI graphs are hidden in the resummation 
and arise via the self consistent Schwinger-Dyson or gap 
equation of the 2PI (see Eq. or 2PPI formalism 

(see Eq. I|3.12l) ). respectively. 

In Fig. ^] we present such a generic 1PI but 2PPR 
graph in the two-loop 2PPI approximation. In Fig. [T51 
we display a similar graph in the two-loop 2PI approach. 
As the 2PPI formalism resums all local contributions to 
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the propagator no ladder diagrams are introduced via 
resummation. In the 2PI formalism in addition nonlocal 
insertions are taken into account which lead to infinite 
ladder resummations. An example for such a ladder di- 
agram is depicted in Fig. It can be identified in the 
lower part of Fig. EH 

As ladder diagrams do not fall apart if two lines meet- 
ing at the same point are cut, they are indeed 2PPI and 
thus join in the 2PPI formalism explicitly as higher order 
corrections to the effective action functional T. We have 
shown a three-loop diagram of "ladder-type" in Fig. 




FIG. 12: A generic 1PI but 2PPR graph, which is produced 
via the resummation in the gap equation for M 2 (t) in the 
two-loop approximation in the 2PPI formalism. Thin solid 
lines denote the free propagator while dashed lines denote 
the classical field 4>. The solid dot indicates the external time 
t. 




FIG. 13: A generic 1PI but 2PR graph, which is hidden in the 
resummation via the full Schwinger-Dyson equation 1A.4> for 
the self energy E' 2 ' (t, t'; p). Thin solid lines denote the free 
propagator while dashed lines denote the classical field <j>. The 
solid dots indicate the external times t and t' . 

The BVA and NLO-l/N approximations in the 2PI 
formalism sum an even larger class of diagrams as in 
these approximations already r 2PI [(f), G] contains an infi- 
nite series of vacuum diagrams with all loop orders. This 
infinite series of chain diagrams can be formulated in 
a very compact way within the auxiliary field formal- 
ism 0, • The diagrams have the topology of chains of 
bubble- graphs (see Fig. ll5l for two generic vacuum graphs 
contributing to T 2PI ). Depending on a given approxima- 
tion these diagrams contribute in the 2PPI formalism as 
well. 



* • 

FIG. 14: Example for a ladder diagram. Thin solid lines 
denote the free propagator while dashed lines denote the clas- 
sical field <f>. The solid dots indicate the external times t and 
t'. 




FIG. 15: Vacuum graphs with a topology of closed chains 
which contribute to the 2PI effective action in the NLO-l/N 
approximation; solid lines represent the 2PI propagator G, 
while the dashed lines denotes the classical field <j>. 



APPENDIX B: SOME MORE COMMENTS ON 
THE NUMERICS AND MOMENTUM 
INTEGRATIONS 

In our simulations we have used a momentum cutoff 
of p max = 20 and an equidistant momenum grid with 
Ap = 0.05. Both choices are somewhat generous; we have 
not attempted an optimization with respect to CPU time 
and storage requirements as one would certainly have to 
do for simulations in 3 + 1 dimensions. In this Appendix 
we more closely investigate the cutoff and momentum 
grid dependences. 

(i) The choice of Ap. We have repeated the simulation 
of Fig.[5]for values of Ap between 0.04 and 0.2 while 
leaving p max = 20 fixed. We display in Fig.lTolthe 
time evolution of the classical field 4>(t) and of the 
effective mass M. 2 (t). The numerical results for 
these quantities are seen to converge for Ap < 0.07. 
The curves for Ap = 0.067 and Ap = 0.04 cannot 
be distinguished. While the qualitative behavior 
of (f>(t) does not change even for larger values of 
Ap the late time averages of the mass M. 2 show a 
considerable dependence beyond Ap ~ .07. 

(ii) The cutoff dependence. The momentum cutoff 
is a cutoff for convergent integrals. As one may 
conclude already from the momentum spectra dis- 
played in Fig. ^2 the cutoff can be reduced appre- 
ciably. In Fig.El we show the dependence of M 2 (t) 
on p max . One sees that even for a cutoff as low as 
Pm&x = 5 the deviations are only at the percent 
level and for p max = 5 the results are already sat- 
isfactory. This may change at later times if the 
momentum distributions get broader by rescatter- 
ing. 
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(a) 

FIG. 16: Detailed study of the dependence on the momentum 
evolution of the mean field <j> (b) evolution for the effective mass 
the distance between the grid points Ap; the solid line represents 
line Ap = 0.1 and the dotted line Ap = 0.2. 
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(b) 

grid for the simulation with the parameters from Fig. |U] (a) 
M 2 ; the momentum cutoff is fixed at p ma x = 20 and we vary 
Ap = 0.04, the long dashed line Ap = 0.067, the short dashed 




tant momentum grid. Indeed it is more economical to 
choose Ap small for small momenta and to let it increase 
for larger ones, as was done, .e.g., in previous computa- 



FIG. 17: Dependence of the time evolution of M 2 {t) on the 
momentum cutoff p max for the simulation with the parameters 
from Fig. [5] in the time range t G [25,50]; in the inset the 
whole time range is shown; the momentum distance is fixed 
at Ap = 0.05 while p max varies between 5 and 20; the solid 
line represents the simulation for p max = 20, the long dashed 
line Pmax = 10, the short dashed line p ma x = 7 and the dotted 
line pmax = 5. 



(iii) The time grid. For our simulations we have cho- 
sen At = 0.001, except for the simulation in 
Pig. IU where At = 0.005. We compare the results 
for the simulation in Fig. obtained with At — 
0.0005,0.001 and 0.005. The results for the first 
two values agree very well; those for At = 0.005 
start to differ at late times. This means that a 
choice At = 0.001 is appropriate. In the case of 
Fig. 8 the variations with time are much slower, so 
that the choice At = 0.005 is sufficient. 

We would finally like to point out that it is in no way 
inherent in our numerical approach to use an equidis- 




FIG. 18: Dependence of the time evolution of M 2 (t) on the 
time At for the simulation with the parameters from Fig. |S]in 
the time range t £ [25, 50]; in the inset the whole time range 
is shown; the momentum distance is fixed at Ap = 0.05 and 
the momentum cutoff p ma x = 10; the solid line represents the 
simulation for At = 0.0005, the long dashed line At = 0.001 
and the short dashed line At = 0.005. 



tions of our group. In one space dimension the choice 
of equidistant momenta turns momentum conservation 
in a trivial algebra of indices. In three space dimensions 
one may use the 0(3) invariance of the mode functions 
as functions of p. Then, due to angular integrations, 
the momentum integrals become convolutions of mode 
functions with phase space functions and an equidistant 
momentum grid does not lead to any major simplifica- 
tion. 
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